Chaos suppression in the parametrically driven Lorenz system 
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We predict theoretically and verify experimentally the suppression of chaos in the Lorenz system 
driven by a high-frequency periodic or stochastic parametric force. We derive the theoretical criteria 
for chaos suppression and verify that they are in a good agreement with the results of numerical 
simulations and the experimental data obtained for an analog electronic circuit. 
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I. INTRODUCTION 

Control of the chaotic dynamics of complex nonlinear 
systems is one of the most important and rapidly develop- 
ing topics in applied nonlinear science. In particular, the 
concept of chaos control, first introduced in Ref. 0, has 
attracted a great deal of attention over the last decade. 
Among different methods for controlling chaos that have 
been suggested by now, the so-called non-feedback con- 
trol is attractive because of its simplicity: no measure- 
ments or extra sensors are required. The idea of this 
method is to change the complex behavior of a nonlinear 
stochastic system by applying a properly chosen external 
force. It is especially advantageous for ultra-fast pro- 
cesses, e.g., at the molecular or atomic level, where no 
possibility to measure the system variables exists. 

Many of the suggested non-feedback methods employ 
the external forces acting at the system frequencies [2|, 
including parametric perturbations with the frequency 
that is in resonance with the main driving force. In par- 
ticular, changing the phase and frequency of a parameter 
perturbation in a bistable mechanical device was shown 
to either decrease or increase the threshold of chaos y|- 
Similarly, unstable periodic orbits are known to be stabi- 
lized by a low-frequency modulation of a system control 
parameter Q, when the control frequency is much lower 
than the system frequency. 

A possibility to change significantly the system dynam- 
ics by applying a high- (rather than low- ) frequency force 
is known for almost a century. As a textbook example, 
we mention the familiar stabilization of a reverse pendu- 
lum (known as the Kapitza pendulum) by rapid vertical 
oscillations of its pivot § . This discovery triggered the 
development of vibrational mechanics [|| where the gen- 
eral analysis of nonlinear dynamics in the presence of 
rapidly varying forces is based on the Krylov-Bogoljubov 
averaging method .7]. In the control theory, the high- 
frequency forces and parameter modulations are usually 
used for the vibrational control of non-chaotic nonlin- 
ear systems However, as was shown for the exam- 
ple of the Duffing oscillator [|| Il3> chaos suppression 
can also be achieved by applying a high-frequency para- 
metric force. Later, chaos suppression in the Belousov- 



Zhabotinsky reaction was demonstrated numerically by 
adding white noise , and the effect of random para- 
metric noise on a periodically driven damped nonlinear 
oscillator was studied by the Melnikov analysis [T3| . 

In this paper we apply the concepts of the nonres- 
onant non-feedback control to the Lorenz system and 
demonstrate analytically, numerically, and also experi- 
mentally that the suppression of the chaotic dynamics 
can be achieved by applying a high-frequency parametric 
or random parametric force. The Lorenz system, found 
in 1963, is known to produce the canonical chaotic at- 
tractor in a simple three-dimensional autonomous sys- 
tem [l3l ll l| . and it can be applied to describe many 
interesting nonlinear systems, ranging from thermal con- 
vection [15j to lasers dynamics jlq . 

The paper is organized as follows. In Sec. II we present 
our model and outline the theoretical method and results. 
By applying the averaging method, we derive the effec- 
tive Lorenz equations with the renormalized control pa- 
rameter and obtain the conditions for chaos suppression. 
In Sec. Ill, we demonstrate the suppression of chaos by 
means of direct numerical simulations and also present 
the experimental data obtained for an analog electronic 
circuit. Finally, Sec. IV concludes the paper. 



II. THEORETICAL APPROACH 
A. Model 

We consider the familiar Lorenz system [l3| driven by 
a parametric force. In the dimensionless variables, the 
nonlinear dynamics is governed by the equations 

x = cr(y- x), 
y = r[l + f{t)]x -y-xz, 



(1) 



xy 



bz, 



where the dots stand for the derivatives in time, a, r, and 
b are the parameters of the Lorenz model, and the func- 
tion f(t) describes a parametric force. For definitcness 
but without restrictions of generality, we select the stan- 
dard set of the parameter values, a = 10 and b — 8/3, 
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whereas the parameter r is assumed to vary. It is well 
known LLi| that, in the absence of the parametric driv- 
ing force f(t), the original Lorenz equations demonstrate 
different dynamical regimes on variation of the control 
parameter r, which are associated with the existence and 
stability of several equilibrium states. In brief, the sys- 
tem dynamics can be characterized by three regimes: 

• the G\ regime, for r < 1: there exists the only 
stable fixed point at the origin, {x, y, z) = (0,0,0); 

• the G2 regime, for 1 < r < 24.74: the origin 
becomes unsta ble, two ne w fixed p oints appear, 
(x, y, z) = (±y/b(r-l), ±y/b(r - Tj/r, r - 1); 

• the G3 regime, for r > 24.74: no stable fixed points 
exist, chaotic dynamics occur with a strange attrac- 
tor. 

We consider the special case of the general model Q , 
assuming that the characteristic frequency of the para- 
metric force f(t) is much larger than the characteristic 
frequency of the unforced Lorenz system, where the fre- 
quency (which is the mean-time derivative of the phase) 
of the Lorenz system can be defined as 17] 



oj 



lim 

T^oo 



2irN(T) 
T : 



(2) 



where N(T) is the number of turns performed in T. 



B. Periodic driving force 



First, we consider the case of a periodic driving force, 



f(t) = k cos(u>t). 



(3) 



Assuming that the frequency of parametric modulations 
is large (uj 3> ujq), we apply an asymptotic analytical 
method 0, |(| based on a separation of different time 
scales, and derive the effective equations that describe 
the slowly varying dynamics. To do this, we decompose 
every variable into a sum of slowly and rapidly varying 
parts, i.e., 



x = X + f , y = Y + r], z = Z + C, 



(4) 



where the functions £(t), T)(t), and ((t) describe fast os- 
cillations around the slowly varying envelope functions 
X(t), Y(t), and Z(t), respectively. The rapidly oscil- 
lating corrections are assumed to be small in compari- 
son with the slowly varying parts, and their mean values 
during an oscillation period vanish, i.e. < x >= X, 
< y >= Y, and < z >= Z. Substituting the Eq. (gj 
into Eq. (Q with the force J3J), and averaging over the 
oscillation period, we obtain the following equation, 



X = a(Y -X), 

Y = rX - Y - XZ + rk < £cos(ujt) >, 



(5) 



where the terms < > and < £77 > are neglected 
because, for large uj, they are all of the higher orders 
in kuj^ 1 . Using Eq. J^J and keeping only the terms 
not smaller than those of the order of kuj -1 and k, 
respectively, we find the equations, £ = arj and 7) = 
rkX cos(wt). Regarding the function X as constant dur- 
ing the period of the function f(t), we obtain the solu- 
tion £ = — ar(k/ui 2 )X cos(wrj), and thus < £cos(a;£) >= 
— ar(k/2uj 2 )X . Therefore, the averaged equations JSJ are 



X = a{Y - X), 

Y = r cH X -Y-XZ, 

Z = XY - bZ, 



(6) 



where 



r e g = r(l - rKu), 



2uj 2 



> 0. 



(7) 



As a result, the averaged dynamics of the Lorenz system 
in the presence of a rapidly varying parametric force is 
described by the effective renormalized Lorenz equations 
© with the effective control parameter r c fj, so that all 
dynamical regimes and the route to chaos discussed ear- 
lier can be applicable directly to Eq. ©, assuming the 
effect of the renormalization. 

Therefore, in terms of the averaged system, for r e g < 1, 
i.e. under the condition 



K u >{r- l)/r 2 



(8) 



the fixed point at the origin remains stable and, therefore, 
in terms of the original Lorenz system, the chaotic dy- 
namics should be suppressed. Next, for 1 < r c s < 24.74, 
i.e. under the condition 



(r - 24.74)/r 2 < < (r - l)/r 2 



(9) 



the Lorenz chaos is also suppressed and the stable saddle- 
focus points appears: X = Y = ±^b(r c g — 1), Z = 
r c s ~ 1- 

Relations JHJ) and JHJ that follow from the averaged 
equations, define the conditions for the chaos suppres- 
sion; they can be expressed as relations between the am- 
plitude and frequency of the rapidly varying oscillations, 



k > u^2{r- l)/ar 2 



w v /2 (>-24.7)/err 2 < k < u^2(r - l)ar 2 



(10) 



(11) 



Z=IF- bZ, 



The dependencies K u ~ r (see Eqs. |JSJ| and J??}) and 
k ~ uj (see Eqs. (|10|) and (1111) ) are the key characteristics 
of the chaos suppression effect that we verify numerically 
and compare with the experimental data obtained with 
an analog electronic circuit (see Figs. 4 and 5 below). 
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C. Random driving force 

Now we turn to the case of a random force, and treat 
the function f(t) in Eq. Q as random by formally writ- 
ing f(t) = e(t), where e(t) describes a bandwidth-limited 
noise with a power spectral density 

fp(u)), when u>\ < \u>\ < u>2, 
(12) 
0, when \w\ < uj\, \ui\ > 

and the zero mean value, < e(t) >= 0. In order to apply 
the analytical method discussed above, the noise e(t) is 
assumed to be composed of high-frequency components 
only, i.e. ui S> wq, where 



u)P e (id)du> 



P e (u))du> 



is the characteristic frequency of the bandwidth-limited 
parametric fluctuations. Again, decomposing the system 
variables into sums of slow and rapidly varying functions 
and averaging over the period To, we obtain 



X 



a{Y-X), 

rX-Y-XZ + rk< £(t)e(t) >, 
Z = XY - bZ. 



Y 



(13) 



where the value of T is much smaller than the charac- 
teristic time scale of the Lorenz system oscillations but 
sufficiently larger than the time scale of fluctuations, i.e. 
2-k/uj < T "C 2tt/uj q . Then, the equations for the rapid 
parts can be reduced to the equations, £ = crr](t) and 
fj = re(t)X. Thus, the key averaged quantity becomes 
< £(t)e(t) >=< £(*)£(*) > /(o-rX). On the other hand, 
for the stationary random process £(t), the following ex- 
pressions are valid, 



< mm >=-< i 2 (t) >-- 



uu 2 P^(ui)du) , 



where Pf(u>) = \H (lS)\ 2 P t (Lo) is the power spectral den- 
sity of £(t) and H(ui) — —crrX/u) 2 is the frequency re- 
sponse function for the equation £ = arXe(t). As a re- 
sult, we obtain the averaged equations in the form ifSjl 
where this time [cf. Eq. {7J|] 



r e g = r(l - rK e ), K e = 2a 



p{w) 



duj. 



(14) 



The quantity K e is related to the intensity of the para- 
metric fluctuations e(t) and it is always positive. Since 
the averaged equations have the same form as the original 
unforced Lorenz equations but with the effective control 
parameter r g instead of r, there appears the parame- 
ter region where the chaotic dynamics is suppressed. In 
particular, for r e g < 1, i.e., under the condition 



K e > (r - l)/r 2 



(15) 
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FIG. 1: (a) Numerical simulation of chaos suppression shown 
for the evolution of the function x(t). The high-frequency 
parametric force is turned on at the point t = 20, and it 
moves the system to a new fixed point < x(t) >= 6.8, for 
k = 3.5, or it stabilizes the origin x(t) = 0, for k — 6.5. (b,c) 
Phase plane representation of the system dynamics for (b) 
k — 3.5 and (c) k = 6.5. 
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FIG. 2: Phase plane representation of the system dynamics 
after the suppression of chaos, at different values of the para- 
metric force amplitude (marked in the figure). 



the fixed point at the origin remains stable in the pres- 
ence of rapidly varying fluctuations. When 1 < r c g < 
24.74, i.e. for 



(r - 24.74)/r 2 < K e < (r - l)/r 2 , 



(16) 



two new saddle-focus fixed points appear, with the coor- 
dinates X = Y = ±-y/6(r B ff - 1), Z = r e fj — 1. 

When the random function e(t) describes a bandwidth- 
limited white noise, P e (u>) — So, the renormalization fac- 
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FIG. 3: Schematic of the analog electric circuit representing 
the Lorenz oscillator in the rescaled variables, including the 
parametric driving force and fluctuations. 

tor K e can be written as 

LU1LO2 LU1LO2 
and the effective control parameter is simplified to be 

r eff = r (l - — < e 2 (t) >) , (18) 

so that the chaos suppression effect is clearly proportional 
to the noise intensity. Moreover, in the limit u>i,u)2 — > u>, 
we recover the result K e = ak 2 /(2a; 2 ) obtained earlier 
for the periodic parametric force [see Eq. Q]. 



III. NUMERICAL SIMULATIONS AND 
EXPERIMENTAL RESULTS 

To verify out theory, first we perform direct numeri- 
cal simulations of the full model (JTJ . Figures ^ an( l [21 
show the results for the system temporal evolution ob- 
tained numerically by solving Eqs. and © with the 
parameters: (<j,b,r) — (10,8/3,28). Since the mean fre- 
quency of the unforced oscillations of the Lorenz system 
is found to be u>o = 8.24 [see Eq. (j2J)], for the above pa- 
rameters the frequency of the rapid parametric force is 
chosen as to = 70. As shown in the examples presented 
in Fig. the fixed point at the origin x — is stabilized 
for k — 6.5 whereas a new stationary point Xq = 6.8 
appears for k = 3.5. We can see from Fig. |2 that the sys- 
tem dynamics in the phase space after the suppression of 
chaos is reduced to the coexisting limit cycles (marked in 
solid and dotted) whose averaged behavior is described 
by Eq. © , shrinking to the origin as the parametric force 
amplitude k is increased. 
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FIG. 4: Experimental vs. theoretical results shown as the 
stability regions in the parameter plane [K^/K^, r). Solid: 
theoretical result [Eq. JHJ and Eq. ll5ll l for stabilization of the 
fixed point at the origin. Dashed: the theoretical results @ 
and 11 1 (it for the transition to a new oscillating state. Squares 
and diamonds: experimental results for the periodic driving 
force with the frequency to = fKHz. Asterisks: the critical 
values K e for the chaos suppression by a bandwidth-limited 
white noise with U\ = 700Hz and L02 = f300Hz. 



Next, we study experimentally the chaos suppression 
in the Lorenz system in the framework of an analog elec- 
tronic circuit implementation of the Lorenz system (see, 
e.g., Refs. PJI [jjj ). Our circuit uses three Op- Amps 
(TL084) as the building blocks for the operations of sum, 
difference, and integration, and three analog multipliers 
(AD633) for the operations of product. 

As the first step, we rescale the variables x, y, and z 
in order to fit within the dynamical range of the source 
(-15V, 15V), and such that the circuit operates in the 
frequency range of a few kHz. Applying the transforma- 
tions u = x/10, v = y/10, w = z/10, and r = t/A, we 
obtain the rescaled Lorenz system, where u>' = Au>, so 
that Lot — uj't (t is measured in seconds). The rescaled 
equations describe the temporal evolution of the corre- 
sponding output voltages of Op- Amps in the analog elec- 
tronic circuit, as shown schematically in Fig. [3J 

In the electric circuit shown in Fig. |3 the analog mul- 
tipliers (AD633) have a noticeable offset at the output 
that may alter the dynamical behavior of the system; this 
effect has been compensated by using a compensation ar- 
ray. The parameters of the Lorenz model can be defined 
through the values of the resistors, a = Rq/R\ = Rq/R2, 
b = Rq/Rs, r — R 6 /R 3 , and the values of the re- 
sistors used are: R% = R2 = lOOfcO, R3 = 35.7kil 
(for r = 28), i? 4 = R 7 = lOfcfi, R 5 = bkQ, and 
R$ = 374fcf2, whereas the values of the capacitors are: 
Ci = C 2 = C 3 = C = WOnF. The tolerance of the 
resistors and capacitors is less than 1%. For these val- 
ues, the parameters take the following values: a = 10, 
b = 2.67, and A = 1/RqC = 10(s" 1 ). The control pa- 
rameter r of the Lorenz system and the force amplitude 
k can be adjusted by varying the resistor R3 through 
the relation r = 10 6 57/i?3 and changing the amplifica- 
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FIG. 5: Experimental vs. theoretical results for the critical 
amplitude of the parametric force vs. frequency for r = 28. 
Analytical results (I10B and 11 U are plotted as solid and 
dashed lines, respectively. Squares and diamonds show the 
experimental results for the stabilization of the origin and 
the suppression of chaos into a fixed point, respectively. 

tion coefficient of the gain-block element, respectively, as 
shown in Fig. |31 In the electric circuit, the frequency 
of the unforced oscillations for r = 28 is found to be 
lu' q = Aqjq = 13.1Hz= 82.4 rad/s. Therefore, we apply a 
parametric force with the frequency larger than 65Hz. 

Figure 01 shows the threshold curves described by 
Eqs. ijSJl and JjJJl, as well as Eqs. 1(15)1 and that define 
three regions with different nonlinear dynamics, as com- 
pared with the numerical and experimental data. The 
domains G\ and Gi indicate the regions for stabilizing 
the fixed point at the origin and the transition to a new 
fixed point, respectively, whereas the domain G3 is the 
region of the chaotic dynamics. The critical values of 



K u were obtained experimentally by changing the am- 
plitude k of the parametric driving force at the fixed fre- 
quency oj — 1kHz, the corresponding experimental data 
are presented in Fig. 01 by squares and diamonds. In the 
experiment with parametric fluctuations (the data are 
marked by asterisks), bandwidth- limited white noise with 
the limits uj\ = 700Hz and ll>2 = 1300Hz was applied, and 
its strength was adjusted by gain-block element (Fig.[3J|. 

Figure [3] shows the dependence of the critical ampli- 
tude of the parametric force for the chaos suppression as 
a function of the frequency. Solid and dashed lines rep- 
resent the theoretical results described by Eq. (|10|) and 
Eq. (|llf> . respectively, whereas the experimental data are 
plotted as squares and diamonds. As shown in Fig.0]and 
Fig. the experimental results are in a good agreement 
with the theoretically calculated threshold dependencies. 



IV. CONCLUSIONS 

We have studied analytically, numerically, and exper- 
imentally the suppression of chaos in the Lorenz system 
driven by a rapidly oscillating periodic or random para- 
metric force. We have derived the theoretical criteria for 
chaos suppression which indicate that, for a fixed value 
of the control parameter r, the critical amplitude of the 
force required for the suppression of chaos is proportional 
to its frequency. The theoretical criteria for chaos sup- 
pression have been found to agree well with both the re- 
sults of numerical simulations and the experimental data 
obtained for an analog electronic circuit. 
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